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Abstract 

A search for the Standard Model Higgs boson in the decay channel H — > ZZ^*^ — £~^£~£~^i~, where 
i = e, /U, is presented. Proton-proton collision data at ^/s = 7 TeV recorded with the ATLAS detector 
and corresponding to an average integrated luminosity of 2.1 fb~^ are compared to the Standard Model 
expectations. Upper limits on the production cross section of a Standard Model Higgs boson with a mass 
between 110 and 600 GeV are derived. The observed (expected) 95% confidence level upper limit on the 
production cross section for a Higgs boson with a mass of 194 GcV, the region with the best expected 
sensitivity for this search, is 0.99 (1-01) times the Standard Model prediction. The Standard Model Higgs 
boson is excluded at 95% confidence level in the mass ranges 191 — 197, 199 — 200 and 214 — 224 GeV. 
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1. Introduction 



The search for the Standard Model (SM) Higgs 
boson [H-Q is a major goal of the Large Hadron 
Collider (LHC) programme. Direct searches at the 



CERN LEP 



collider led to a lower limit on the 



Higgs boson mass, rrin, of 114.4 GeV at 95% confi- 
dence level (CL) 4]. The searches at the Fermilab 
Tevatron pp collider have excluded at 95% CL the 
region 156 GeV < mn < 177 GeV [sj. Results 
from the 2010 LHC run extended the search in the 
region 200 GeV < mn < 600 GeV by excluding a 
Higgs boson with cross section larger than 5 — 20 
times the SM prediction [6, 7]. 

This letter presents a search for the SM Higgs 
boson in the mass range from 110 to 600 GeV 
in the channel H ZZ^*^ £+£-£'+£'-, where 
£,£ — e, II. Three distinct final states, fifJ-fifJ- (4/^), 
ee/ifi (2e2^), and eeee (4e), are selected. The 
largest background to this search comes from con- 
tinuum ZZ^*^ production. For mn < 180 GeV, 
contributions from Z -|- jets and tt processes, where 
the additional charged leptons arise either from 
semi-leptonic decays of heavy flavour or from light 
flavour jets misidentified as leptons, are impor- 
tant. The pp collision data were recorded with the 
ATLAS detector at the LHC at ^/s = 7 TeV and 



correspond to an average integrated luminosity of 
2.1 fb-i 8]. 



2. The ATLAS Detector 

The ATLAS detector [9] is a multi-purpose parti- 
cle physics apparatus with forward-backward sym- 
metric cylindrical geometrj0. The inner tracking 
detector (ID) consists of a silicon pixel detector, 
a silicon microstrip detector, and a transition ra- 
diation tracker. The ID is surrounded by a thin 
superconducting solenoid providing a 2 T mag- 
netic field. A high-granularity lead-liquid argon 
(LAr) sampling calorimeter measures the energy 
and the position of electromagnetic showers. An 
iron-scintillator tile calorimeter provides hadronic 
coverage in the central rapidity range. The end- 
cap and forward rapidity regions are instrumented 
with LAr calorimetry for both electromagnetic and 
hadronic measurements. The muon spectrometer 



^ ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point. The z-axis is along 
the beam pipe, the x-axis points to the centre of the LHC 
ring and the y-axis points upward. Cylindrical coordinates 
(r,</>) are used in the transverse plane, <j> being the azimuthal 
angle around the beam pipe. The pseudorapidity rj is defined 
as T) = — ln[tan(6/2)] where 6 is the polar angle. 
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(MS) surrounds the calorimeters and consists of 
three large superconducting toroids, each with eight 
coils, a system of precision tracking chambers, and 
detectors for triggering. A three-level trigger sys- 
tem selects events to be recorded for offline analysis. 



3. Data and Simulation Samples 

The accumulated data are subjected to quality 
requirements ensuring that the relevant detector 
components were operating normally. The result- 
ing average integrated luminosity of 2.1 fb^^ cor- 
responds to 2.28 fb-\ 1.96 fb~i and 1.98 Ad^^ for 
the 4/i, 2e2/i and 4e final states, respectively. 

The H ZZ'-*^ M signal is modelled us- 
ing the POWHEG Monte Carlo (MC) event genera- 
tor Hi, which calculates separately the gluon 
and vector-boson fusion production mechanisms 
with matrix elements up to next-to-leading order 
(NLO). The Higgs boson transverse momentum, 
Pt, spectrum is reweighted to the calculation of 
Ref. |12| . providing QCD corrections up to next- 
to-leading order and QCD soft-gluon resummations 
up next-to-next-to- leading log (NNLL). powheg 
is interfaced to pythia [l^l for showering and 
hadronization, which in turn is interfaced to PHO- 
TOS 



14l | for QED radiative corrections in the final 
state and to tauola 1^ 1^ ioi the simulation of 
r decays. 

The cross sections for Higgs boson production, 
the corresponding branching fractions, as well as 
their uncertainties [l7| , are derived to next-to-next- 
to-leading order (NNLO) in QCD for the gluon fu- 
sion 1814231 and vector boson fusion [24] processes. 
In addition, QCD soft-gluon resummations u p t o 
NNLL are available for the gluon fusion process [25| , 
while the NLO electroweak (EW) corrections are 
applied to both the gluon fusion [H, and vector 
boson fusion H, 2^ processes. The Higgs boson de- 
cay branching ratio to the four-lepton final state is 
predicted by prophecy4f j3C), iil] , which includes 
the complete NLO QCD-I-EW corrections, interfer- 
ence effects between identical final state fermions 
and leading two-loop heavy Higgs boson corrections 
to the four-fermion width. Table [T] gives the pro- 
duction cross sections for the iJ — > 4^ for several 
Higgs boson masses. 

The ZZ^*") background is generated using 
PYTHIA, taking into account Z — 7 interference. For 
the inclusive total cross section and the shape of the 
mzz(') spectrum, the mcfm [s^, HI] prediction is 



Table 1: Higgs boson production cross sections for both 
gluon and vector-boson fusion processes in pp collisions at 
y/s = 7TeV. The cross sections include the branching ratio 
of -ff — > A£, with £ = e, fi. The errors arc the total theoretical 
systematic uncertainty. 
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used, which includes both quark-antiquark annihi- 
lation at QCD NLO and gluon fusion. The inclu- 
sive Z boson production, Z -I- jets, is modelled using 
ALPGEN [sl] and is divided into Z+light fiavour 
jets and Zhh] overlaps between the two samples 
are removed. Specifically, bb pairs with separation 
A_R = Y^A02 + A772 > 0.4 between the 6-jets are 
taken from the matrix-element calculation, whereas 
for Ai? < 0.4 the parton-shower jets are taken. 
PYTHIA is also used as a cross-check of the ALPGEN 
results. In this search the Z + jets production is 
normalized from the data, but for comparisons the 
QCD NNLO FEWZ [si, and the mcfm [H, [H 
cross section calculations are used for the inclusive 
Z boson and the Zbb production, respectively. The 
tt background is modelled using MC@NLO [37| and 
is normalized to the approximately NNLO cross sec- 
tion calculated using hathor [38]. Both ALPGEN 
and MC@NLO are interfaced to herwig [39] for par- 
ton shower hadronization and to jimmy [40] for the 
underlying event simulations. 

All generated events undergo a full detector sim- 
ulation performed using GEANT4 



411,142]. 



The number of pp interactions in the same bunch 
crossing (pileup) is included in the simulation. The 
MC samples are reweighted to reproduce the ob- 
served distribution in the data. 



4. Physics Object Identification and Event 
Selection 

The data considered in this analysis were se- 
lected using single-lepton triggers. For electrons 
the threshold on the transverse energy. 
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20 — 22 GeV depending on the LHC instantaneous 
luminosity and for muons the threshold on px was 
18 GeV. Both triggers are more than 99.5% efficient 
for events passing the offiine selection described be- 
low. 

Electron candidates consist of clusters of energy 
deposited in the electromagnetic calorimeter asso- 
ciated to ID tracks. The electrons must satisfy 
the "medium" electron criteria |43j], which require 
the shower profiles to be consistent with those ex- 
pected for electromagnetic showers and a well re- 
constructed ID track pointing to the corresponding 
cluster. The electron transverse momentum is com- 
puted from the cluster energy and the track direc- 
tion. 

Muon candidates are reconstructed by matching 
ID tracks with either full or partial tracks in the 
MS [i^. For the former case, the two independent 
momentum measurements are combined, whereas 
for the latter case the momentum is measured us- 
ing the ID information only, with the MS providing 
muon identification. To reject cosmic rays, tracks 
are required to be consistent with having origi- 
nated from the primary vertex, defined as the recon- 
structed vertex with the highest ^ p\ of associated 
tracks. 

Leptons from Higgs boson decays are expected to 
be isolated and to originate from a common vertex. 
Track and calorimeter isolation as well as trans- 
verse impact parameter significance requirements 
are therefore applied to further reduce the Z -I- jets 
and tt contributions. The sum of of tracks within 
Ai? < 0.2 of the lepton divided by the lepton is 
required to be less than 0.15, while the sum i?T of 
the calorimeter cells within Ai? < 0.2 around the 
lepton divided by the lepton is required to be less 
than 0.3. In the case of electrons, the calorimeter 
cells corresponding to the electromagnetic shower 
are subtracted. The transverse impact parameter 
significance, defined as the transverse impact pa- 
rameter of the lepton with respect to the primary 
vertex divided by its uncertainty, for the two low- 
est px leptons of the quadruplet in events with 
mu < 190 GeV is required to be less than 3.5 and 6 
for muons and electrons respectively. The selection 
efficiency of the isolation and impact parameter re- 
quirements has been studied using data both for iso- 
lated leptons, with Z ^ It decays and non-isolated 
leptons from semi-leptonic 6 and c— quark decays in 
a heavy-fiavour enriched dijet sample. Good agree- 
ment is observed between data and simulation. 

Higgs boson candidates are searched by select- 



ing two same-fiavour, opposite-sign isolated lep- 
ton pairs in an event. Each lepton must satisfy 
Pt > 7 GeV and be measured in the pseudorapid- 
ity range |?7| < 2.47 for electrons and < 2.5 for 
muons. The electron px threshold is increased to 
15 GeV in the transition region between the bar- 
rel and end-cap calorimeters (1.37 < |?7| < 1.52). 
At least two leptons must have pt > 20 GeV. 
The leptons are required to be well separated from 
each other with AR > 0.1. The invariant mass 
of the lepton pair closest to the nominal Z bo- 
son mass (mz) is denoted by mi2 and it is re- 
quired that \mz — mi2\ < 15 GeV. The invari- 
ant mass of the remaining lepton pair, 777,34, is re- 
quired to be lower than 115 GeV and greater than a 
threshold depending on the reconstructed four lep- 
ton mass, 7T74^, as summarized in Tabled The fi- 
nal discriminating variable is mu, where the Higgs 
boson production would appear as a clustering of 
events. The width of the reconstructed Higgs bo- 
son mass distribution is dominated by experimen- 
tal resolution at low mn values, with a full-width 
at half-maximum (FWHM) which varies according 
to decay mode and is between 4.5 {4fi) and 6.5 (4e) 
GeV for mn = 130 GeV. At high rriH the recon- 
structed width is dominated by the natural width 
of the Higgs boson with a FWHM of approximately 
35 GeV at mn = 400 GeV. 

5. Background Estimation 

The dominant ZZ^*^ background is estimated us- 
ing MC simulation. Generated events are required 
to pass the complete analysis selection and the final 
yield is normalized to the integrated luminosity. 

The tt background is also estimated using MC 
simulation. Comparison of data to MC predictions, 
in a control sample of events with opposite sign 
electron-muon pairs consistent with the Z boson 
mass and with one or two additional charged lep- 
tons, are used to verify that the ti background is 
small with respect to the dominant ZZ*^*-* process 
and in agreement with expectation. 

The Z + jets background is normalized using 
data. The control sample is formed by selecting 
events with a pair of same-fiavour, opposite-sign 
isolated leptons consistent with the Z boson mass, 
\mz — 7T7i2| < 15 GeV, and a second same-flavour, 
opposite-sign lepton pair where only kinematic, but 
no isolation or impact parameter, requirements are 
applied. At this stage, the dominant background 
source depends on the flavour of the second lepton 
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Table 2: Thresholds applied to 71134 for reference values of m,4i (see text). For other m^i values, the selection requirement is 
obtained via linear interpolation. 
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Table 3: The expected numbers of background events, with their systematic uncertainty, separated into "Low mass" (m^i < 
180 GeV) and "High mass" {ni4£ > 180 GeV) regions. The expected numbers of signal events for different mu hypotheses 
and the observed numbers of events are also presented. 
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Figure 1: Invariant mass distributions [(a)] mi2 , [(b)] m34, and |(c)| r)i4£ for the selected candidates. The data (dots) arc compared 
to the background expectations from the dominant ZZ^*^ process and the sum of tt, Zbb and Z+light flavour jets processes. 
Error bars represent 68.3% central confidence intervals. 



pair: Z+light flavour jets dominates the final states 
with a second electron pair, while Zbh production 
dominates the final states with a second muon pair 
after the contributions from tt, ZZ^*\ and muons 
from in-flight tt and K decays which correspond 



to 44% of the event yield are subtracted. The ob- 
served background, which is found to be in good 
agreement with expectation, is extrapolated to the 
signal region by means of the MC simulation. 
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6. Systematic Uncertainties 

Uncertainties on lepton reconstruction and iden- 
tification efficiency, and on tlie momentum reso- 
lution and momentum scale are determined using 
samples of W, Z and J/ip decays. The muon effi- 
ciency uncertainty results in an acceptance uncer- 
tainty on the signal and the irreducible background 
which is uniform over the mass range of interest and 
amounts to 1.7% (1.2%) for the 4yU {2e2fi) channel. 
The uncertainty on the electron efficiency results in 
an acceptance uncertainty of 3% (2%) for the 4e 
(2e2/Lt) channel at = 600 GeV reaching 15% 
(6%) at m4e = 110 GeV. 

A conservative theoretical uncertainty of 15% is 
assigned to the ZZ^*^ background contribution 44 1. 
The Z-|-light flavour jets and Zbb backgrounds are 
evaluated using data. A systematic uncertainty be- 
tween 20% and 40% is assigned on their normaliza- 
tion to account for the statistical uncertainty in the 
control sample and the MC-based extrapolation to 
the signal region. The uncertainty on the tt cross 
section is found to be 10% by adding linearly the 
contributions from variations of the renormalization 
and factorization scales to those of the parton dis- 
tribution functions. 

The theoretical uncertainties on the Higgs boson 
production cross section are 15-20% for the gluon 
fusion process and 3-9% for the vector-boson fusion 
17 1, depending on the Higgs boson masfl 



process 



They include uncertainties on the QCD scale and 
on the parton distribution functions 46-4^. An 
additional 2% uncertainty is added to the signal se- 
lection efficiency due to the modelling of the signal 
kinematics. This is evaluated by comparing sig- 
nal samples generated with PYTHIA and the default 
POWHEG samples. 

The overall uncertainty on the total integrated 
luminosity is 3.7% Q. 

7. Results 

The number of events observed in each final state, 
separately for m4i < 180 GeV and > 180 GeV, 



^The limits presented in this study for niH > 200 GeV 
assume cross sections based on on-shell Higgs boson produc- 
tion and decay and use MC generators with an ad-hoc Breit- 
Wigner Higgs line shape. Recently potentially important ef- 
fects related to off-shell Higgs boson production and interfer- 
ence effects between the Higg s boson signal and backgrounds 
have been discussed llTl |45|1 . The inclusion of such effects 
may affect limits at high Higgs masses (rriH > 400 GeV). 



are compared with the expectations for background 
and signal for various mn hypotheses in Table |31 In 
total 27 candidate events are selected by the anal- 
ysis: 12 4fi, 9 2e2/i, and 6 4e events, while in the 
same mass range 24 ± 4 events are expected from 
the background processes. The mi2, 77134, and m^i 
mass spectra are shown in Fig. [TJ The m^g distri- 
bution for the total background and several signal 
hypotheses is compared to the data in Fig. [2l The 
selected events have been examined visually and no 
evidence for reconstruction problems was identified. 

Upper limits are set on the Higgs boson cross sec- 
tion at 95% CL, using the CLg modified frequen- 
tist formalism [50| with the profile likelihood test 
statistic [Hlj. The test statistic is evaluated with 
a maximum likelihood fit of signal and background 
models to the observed m^e distribution. Figure [3] 
shows the expected and observed 95% CL cross sec- 
tion upper limits as a function of mn and Table S] 
summarizes the numerical values for selected mn 
points. The consistency with the background-only 
hypothesis is quantified using the p- value, the prob- 
ability that a background-only experiment fluctu- 
ates more than the observation. The most signifi- 
cant deviation from the background-only hypothe- 
sis is observed for niH = 242 GeV with a p- value of 
4.9%. These results do not account for the so-called 
"look-elsewhere" effect js^]. The SM Higgs boson 
is excluded at 95% CL in the mass ranges 191 — 197, 
199 - 200 and 214 - 224 GeV. 
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Figure 2; distribution of the selected candidates, com- 

pared to the background expectation. Error bars represent 
68.3% central confidence intervals. The signal expectation for 
three mn hypotheses is also shown. 
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Table 4: Median expected and observed 95% CL upper limits 
on the Higgs boson production cross section for several Higgs 
boson masses, divided by the expected SM Higgs boson cross 
section. 
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Figure 3: The expected (dashed) and observed (full line) 
95% CL upper limits on the Higgs boson production cross 
section as a function of the Higgs boson mass, divided by 
the expected SM Higgs boson cross section. The green and 
yellow bands indicate the expected sensitivity with itlcr and 
±2(7 fluctuations, respectively. 



8. Summary 

A search for the Standard Model Higgs boson 
in the decay channel H — >• ZZ^*'> — >• Ai based on 
2.1 fb~i of data recorded by the ATLAS detec- 
tor at y/s = 7 TeV during the 2011 run, has been 
presented. No significant excess of candidates is ob- 
served in the mass range between 110 and 600 GeV 
with respect to the expected SM background. The 
observed (expected) 95% CL upper limit on the 
Higgs boson production cross section, in units of the 
SM cross section, is 0.99 (1.01) for mn = 194 GeV, 
the region with the best expected sensitivity for this 
search. The SM Higgs boson is excluded at 95% 
CL in the mass ranges 191 - 197, 199 - 200 and 
214- 224 GeV. 
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